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q^ Abstract 

We consider the velocity dependence arising from scattering through dark multipole moments, 
and its effects on the consistency of the signals observed by DAMA and CoGeNT with the dark 
matter hypothesis. We focus on the effects of the experimental uncertainties on the fits, and show 
that the two experiments combined favor dark matter scattering with a velocity-dependent cross- 
section over standard velocity and spin-independent scattering. When appropriate uncertainties 

in 

are taken into account, we show that agreement of the two signals with each other and with the 

m 

l/~) results of null experiments can be obtained. 
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I. THE DAMA-COGENT PUZZLE 

Recently, CoGeNT pQ has reported an excess of events at low nuclear recoil. It is unknown 
what the source of the events is, but it is consistent with a light Dark Matter (DM) particle 
in the 7-11 GeV range interacting with their germanium crystals. The region is remarkably 
close to the mass and cross-section preferred by the DAMA experiment's [2J observation of an 
annual modulation pHH]. The DM mass window in the several GeV range is well-motivated 
by several models of DM [Tl,[8] , notably models that solve the baryon-DM coincidence problem 

At the same time, the preferred cross-section and DM masses of the two experiments, while 
being close to each other in parameter space, are naively inconsistent with each other. They 
are also naively inconsistent with the results of other null experiments, the most relevant 
of which are the silicon run of CDMS pll CLE] and XENON10 [12]. The CoGeNT region 
overlaps neither with the DAMA region for channeled scattering off iodine, nor with the 
DAMA region for unchanneled scattering off sodium (T3J [Hj. 

However, we have also begun to learn more recently about the experimental uncertainties 
which must be properly taken into account to make a conclusion about the relevance of 
experimental constraints for excluding the light DM window, as well as for the consistency 
of the DAMA and CoGeNT signals with each other. For example, a new measurement of the 
scintillation light yield efficiency, £ e ff [IE], opens an allowed region for DAMA channeled spin- 
independent scattering off iodine [13]. In addition, by choosing the sodium and germanium 
quenching factors appropriately, the DAMA unchanneled and CoGeNT spin-independent 
scattering regions can be marginally consistent with each other [16J. Lastly, by assuming a 
20% systematic uncertainty in the energy threshold for the constraint from CDMS silicon, 
the region in which the two signals are consistent with each other can be made consistent 
with the null observation of CDMS. 

While obtaining consistency in this way is marginally possible for spin-independent 
scattering, it does stretch the experimental parameters and theoretical uncertainties to their 
limits in order to allow it. The purpose of this paper is two-fold. First, we unpack the results 
of [16] to show more explicitly how the experimental and theoretical uncertainties allow for 
better agreement between the signals of DAMA and CoGeNT with each other and with the 
results of null experiments for standard spin-independent scattering. We find, however, that 
pushing all the experimental systematic uncertainties beyond the edge of their la preferred 
values only allows agreement between DAMA and CoGeNT in the 99% C.L. regions. We 
then demonstrate that with the choice of a different operator to mediate the scattering, 
significantly better agreement can be obtained. Tension with the null results of CDMS-Si 
and XENON can also be alleviated. 

Interestingly, the operators that we consider, the dark anapole and magnetic dipole 



moment operators, 

O a = Xl^llXK (1) 

O d = x^xF^/A, (2) 

are unique in that the contributions from spin-dependent and spin-independent scattering 
can be equal for some elements (sodium in particular). 1 The model that we have in mind 
is a massive dark photon kinetically mixed with the visible photon. That the coupling to 
nuclei in the scattering goes through the SM photon imposes constraints on the coefficients 
of the scattering cross-section which we utilize. 

These operators also have unusual velocity and momentum dependence: 
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with M the mediator mass, A an expansion parameter (associated in some models with 
strong coupling in the DM state), ttin the nucleus mass, m x the DM mass, and q and 
v the momentum transfer and velocity of the incoming WIMP. We use the standard 
notation for the form factors Fi,F 2 for a coupling of a gauge field to N, e.g. On oc 
iA^N [Fi'j^ + 2m~ cr i J - l/< i t ') N when A^ is spin-1/2. This unusual momentum and velocity 
dependence has been noted before in other contexts [HJ [TTH23] . though in most of these 
cases only some of the terms in the full expression are considered (but see [8]). We find, 
by contrast, that both terms arising from the magnetic and electric form factors can be 
important and give rise to significantly modified spectra. 

In this paper we show that non-standard velocity and momentum dependence can, 
depending on how they enter into scattering cross-section, reconcile the DAMA and CoGeNT 
regions. The dark magnetic dipole moment interaction in particular has the right structure 
to give agreement between the two experiments, consistent with null results of other direct 
detection experiments. The dark anapole interaction on the other hand does not bring the 
two experimental regions together, and its main benefit is to alleviate tension between DAMA 
and the null results. The magnitude of the shifts in the preferred DAMA and CoGeNT 
regions, and whether this leads to better agreement, is a detailed numerical question. 
This can however be understood qualitatively as follows. CoGeNT records slightly lower 
momentum transfer than DAMA, and since these operators are momentum suppressed, this 
causes CoGeNT to shift slightly up relative to DAMA in comparison to the standard spin- 
independent case. More importantly for these operators, however, is the velocity dependence. 



1 The operator which is usually called the anapole couples to the current, O a = Xl^lsX^iyF^, as discussed 
in [T7]. This operator has the same spin structure as Eq. (Tl]), but has an additional q 2 suppression. 



The maximum momentum transfer depositable in the detector is 2[Anv, where jj,n is the DM- 
nucleus reduced mass. Since the typical momentum transfers observed by CoGeNT are lower 
than those observed by DAMA, and the reduced mass of germanium is higher than that of 
sodium, the typical velocities of particles observed by CoGeNT are significantly lower than 
those of DAMA, further suppressing the overall scattering rate in CoGeNT, and further 
shifting the CoGeNT region up relative to DAMA. The relative importance of these effects 
and the amount that it improves the agreement between DAMA and CoGeNT is a question 
we address in detail in this paper. 

The types of models in which the anapole and magnetic moment operators dominate 
the scattering are not hard to construct. For a Majorana particle scattering through a 
vector mediator, for example, O = xi^X^l^ vanishes, and one expects the anapole to 
be the dominant contribution. For Dirac particles, one can explicitly construct models 
where the coupling of the DM particles is purely axial. Likewise, the DM can have a 
large magnetic moment when constituents charged under a dark force are bound into a 
neutral state. Because the rates are velocity and momentum suppressed, the corresponding 
scattering cross-section must be large. This can be accommodated with a light mediator 
which is weakly coupled to Standard Model particles. We discuss a model where the large 
cross-section generates both the observed rate and is consistent with the results of null 
experiments. 

The outline of this paper is as follows. In the next section we lay out the rates for standard 
spin-independent scattering and for the anapole and magnetic dipole operators. We then 
turn to discussing the effect of the experimental uncertainties on standard spin-independent, 
anapole and magnetic dipole scattering cases, and show how, properly accounting for these 
uncertainties, we can bring the two results into better agreement. We focus in particular on 
the sodium quenching factor, xenon prompt photon (Si) to nuclear recoil conversion factor 
£ e ff, the stochasticity of photo-electrons in XENON10, and the systematic uncertainty in 
the CDMS-Si energy threshold. Lastly, we discuss models that generate the observed event 
rates, and conclude. 



II. SCATTERING RATES 

We begin by reviewing the standard scattering rates and then turn to a discussion of the 
anapole and dipole rates. The rate for scattering is 

dR Ar p x f , 3 „._ _ . da 

- N T ^ d 3 vvf(v,v e )-—, (5) 
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and ix m is the reduced mass of the nucleus-dark matter system. We take the velocity 
distribution f(v,v e ) to be a modified Boltzmann distribution: 

f(v,v e ) oc ( e -C*+W«8 - e-^«) 0(vL - (v + v e f). (7) 

The additional term is to allow for a smooth cut-off of the velocity distribution near the 
galactic escape velocity v esc . The Earth's speed relative to the galactic halo is v e = v + 
f or .&cos7cos[a;(£ — 1 )} with v Q = v + 12 km/s, t> or & = 30 km/s, cos 7 = 0.51, t = June 2nd 
and u = 27r/year. We take as a standard case Uq = 220 km/s but allow vq to vary up to 270 
km/s, with higher values decreasing the tension with XENON10. We fix v esc = 500 km/s; 
allowing this to vary changes the DAMA and CoGeNT windows only by a small amount, 
though the constraints from XENON10 become more stringent for larger v esc . 

A standard calculation relates the differential rate for scattering off nuclei to the scattering 

rate off a nucleus <jn, 

da _ m N a N 

dE R ~ 2/i> 2 ' U 

For the standard spin-independent case, this rate is related to a scattering off protons, cr p , 
through 

T _» 2 N [f P Z + fn(A-Z)} 2 
ftn Jp 

where fx n is the DM-nucleon reduced mass and f p and f n are the dark matter couplings to 
the neutron and proton. We set f n — since this is the normalization that will arise most 
naturally in the models we discuss later, 2 and we consider two different choices for the form 
factor F 2 (Er), which can give rise to O(20%) variations in the derived cross-sections. We 
make use of a Helm form factor 

F(E R ) = 3 Mvo) e-M***,* (10) 

with two different choices for r : 

r =((1.2A 1 / 3 ) 2 -5 S 2 ) V2 fm (11) 

with s — 1, and 

r = ((1.23A 1 / 3 - a) 2 + 7 -7t 2 b 2 - 5sA fm, (12) 

with a = 0.6, b = 0.52 and s = 0.9 E 



2 Choosing f p = f n would give slightly different results. In particular, A 2 /Z 2 is approximately 30% larger 
for germanium than for silicon, and thus the CDMS-Si constraints on CoGeNT would be weakened by 
approximately this amount. Similarly, the region of parameter space favored by sodium scattering at 
DAMA would move up about 20% in cross-section relative to CoGeNT, which would help (hurt) agreement 
between the two regions in the case of scattering through the dipole or anapole (standard) coupling. 



We next derive the rate for scattering through the anapole operator, Eq. (pj). The photon 
coupling to nuclei is 

O n = A"N(p) (F 1 (q)(p + p'), + (F t (q) + F 2 (q))2iE flv q") N(p'). (13) 

where Fi(q), i*2(g) are form factors, and the spin tensor T,^ is a generator in the appropriate 
representation of the Lorentz group for spin- J nuclei N. For instance, for spin-1/2 nuclei, 
^iiv = \vw an d f° r spin-0 nuclei S MJy = 0. In (13), the fields N have the standard non- 
relativistic normalization, which for spin-1/2 nuclei differs from the standard relativistic 
normalization by a factor of ^J2ra^. The form factors satisfy -Fi(O) = Z, (-Fi(O) + -^2(0)) = 
27 m^fe^' wnere b^ denotes the nuclear magnetic moment and b n = e/2m p denotes the 
Bohr magneton, since we are already using the more common symbols fiN,Hn for reduced 
masses. 3 In the non-relativistic limit, the nuclear magnetic moment coupling can be written 
6 at j • B. We take the g-dependence of Fi(q) from the Helm form factor, and we neglect the 
g-dependence of -^2(0). Making these substitutions, the resulting matrix element, for a Dirac 
state, is 

lE W = ^ U^F(E R f - i f fecW '■ z 2 F{ERf _ u ,i£J$!X) . 

4^' ' M 4 V \ m 2 x m 2 N v ' 3J m 2 N b 2 n J J 

(14) 
The resultant scattering cross-section, the analogue of Eq. ^ and which should be inserted 
in Eq. [8] to obtain the differential rate, is 

** = An (^Z*F{E R y - g 2 ( {m *t m 2 N)2 Z*F(E R y - 2A^^-)) . (15) 
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When reporting cross-sections for the anapole case, we use a convention closely related to 



Eq. ([15]), taking a = /i 2 /4vrM 4 . 

Similarly, the rate through the magnetic moment operator, Eq. pi), can be computed. 
We find the resultant scattering cross-section is 

°n = 4tL (4q 2 v 2 Z*F(E R f - g 4 ( ( -*— + ±) Z 2 F(E R f -2A* J+l ^ 
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'(16) 
When reporting cross-sections, we use the convention a = 4/i 2 /7rM 4 . 

In the appendix we offer analytic expressions for the velocity integrals in Eq. (J5l) necessary 
for computing the total rates in both the standard case and in the case of v 2 dependence in 
the rate. We next discuss our results using these expressions for the anapole operator with 
experimental uncertainties folded in. 



3 We are assuming here that the coupling to the nucleus goes through the photon. For a more general 
coupling through a dark force only, the magnetic moment and charge can be allowed to float, shifting our 
results. 



III. RESULTS 

Neither DAMA nor CoGeNT measures the total nuclear recoil energy Er. Instead, both 
experiments measure the "electron equivalent" energy, E ee , which is energy deposited into 
electrons measured as scintillation or ionization. To extract the nuclear recoil energy, and 
hence the energy spectrum of the dark matter recoils, one must fold in the quenching factor 
Q which relates the two. These quenching factors can be highly uncertain. In germanium, 
relevant for CoGeNT, we take two possible parameterizations relating the nuclear recoil and 
electron-equivalent energies 

, E n 0.1204 

Q Ge = 0.19935A Q \-^j , (17) 

and 

/ p \ 0.1204 

Q Ge = 0.224A2 ( -gL J £-Vo-oo383 keV/£* (18) 

Aq is an amplitude for the quenching factor that is allowed to vary between 0.85 and 1 in 
accordance with the uncertainty in the observed quenching factors. 

For DAMA's sodium iodine, the quenching factors carry larger uncertainties. Though 
DAMA reports a quenching factor Q Na = 0.30 ± 0.01 averaged over 6.5 to 97 keV nuclear 
recoil, other measurements seem to suggest that the measured error on the quenching factor is 
much larger. For example, [25J reports a sodium quenching factor Q^ a — 0.4 ±0.2 over 5-100 
keV. Other experiments report smaller quenching factors with smaller errors: 0.275 ± 0.018 
over 4-252 keV |26J, and 0.25 ± 0.06 at 10 keV [27]. On the other hand, in many cases it 
is anticipated that the quenching factor will rise at the lower recoil energies associated with 
DAMA, as shown in [25] . We will take the standard value QNa — 0.3, but show the effect of 
allowing it to rise as high as 0.45 at 5 keV. 

We first show in Fig. (fTl) the results for the standard case, with the quenching factors 



QNa — 0.3 and Qj = 0.09 and the Germanium quenching factor parametrized by Eq. (17). 



We take Vq = 220 km/s, v esc = 500 km/s and the standard Helm form factor with the 



standard r given by Eq. (11). One sees that poor agreement is obtained between CoGeNT 
and unchanneled DAMA 4 . One can also see that poor agreement between DAMA and the 
null results of CDMS-Si is obtained. For the CDMS-Si curves, we take the constraints derived 
from the set-up of [TUlfTT] . and we derive results from the SIMPLE experiment for 14 kg-days 



Obtaining constraints from XENON10 is more involved. Recently, it was pointed out 
that the effects of stochasticity in the number of Si photo-electrons (PEs) observed can be 
important for constraining DM in the low mass window [21] • In particular, each nuclear 



Wc focus only on unchanneled DAMA in this discussion, following 



recoil energy can be mapped onto the number of prompt S\ PEs expected in each scatter. 
To suppress background contamination, XENON10 cuts out events with fewer than two PEs, 
which for a given £ e g roughly corresponds to a recoil energy threshold E^hr given implicitly 
by Ethr = 0.37/C e s(E R ). However, there can be low-energy recoil events below E t \ a that, 
thanks to Poisson fluctuations and detector resolution, nevertheless make it over the SI cut. 
Since the DM signal is highly peaked near low energies, this results in more DM signal 
events being tipped above threshold by the fluctuations than below, resulting in a tighter 
constraint. This was utilized by XENON100 to purportedly rule out the CoGeNT window 
[3T] . In that analysis, however, a high £ e ff converting expected S\ photo-electrons to nuclear 
recoil energy was employed. A choice of a larger £ e s translates to a lower energy threshold 
and a strong constraint on light DM. We think it likely that the £ e g chosen there is too 
high, and hence the constraints are too strong. Here we focus on XENON10 constraints, 
which are in general more restrictive than those due to XENON100 on account of its higher 
threshold. 

We will follow the procedure outlined in [32] , and make use of the detector resolution and 
efficiency derived in that analysis. In particular, we show in each constraint plot the 90% 
limit based on two different assumptions on £ e ff. Both assumptions take the central values 
measured for C e g in [15] , but the first (stronger) constraint assumes a constant £ e ff below the 
lowest measured energies, and the second (weaker) constraint assumes a linear interpolation 
to zero at vanishing recoil energy. The constraint on a WIMP model is determined by the 
predicted rate of events after convolving with the detector resolution and taking into account 
detector and cut efficiencies: 

dR f , ,_, dR , ,_, . dN , „ . . „ . . , 

/ dE R — r r 1 {E R )-—{E R )e{E obs ). (19) 

J2keV aiS R U-^obs 



dE b s 

Here, r\ and e are detector and cut efficiencies (including the acceptance rate of the "50%" 
acceptance box), taken from [32] and [12] respectively, and dN/dE obs is the detector 
resolution also from [32], which is somewhat broader than a pure Poisson distribution as 
utilized in [33]. We further cut off the low end of the integration over recoil energies at 
E R = 2 keV, though due to the low detection efficiency at such low energies, the constraint 
is not strongly dependent on the exact limit. The constraints derived from these two 
extrapolations of C e g are shown as the edges of the green band in the figures. 

Agreement between DAMA and CoGeNT can be improved by choosing a different set 
of experimental parameters. In the second panel of Fig. (1), we show the results with 



Qno, = 0.45 and Qj = 0.09 and the Germanium quenching factor parametrized by Eq. (18). 
Some improvement can be obtained by shifting the quenching factor amplitude in Eq. (18) 
Aq to 0.85 as shown in Fig. pj, left panel. We also modify the form factor r from Eq. (12), 
and take into account a possible 20% uncertainty in the energy threshold of CDMS-Si [3? . 
This is shown in Fig. d2|, right panel. However, even with this extreme set of parameters, 
only marginal agreement can be found between the two results. 



We next consider velocity dependent cross-sections, and we look at the anapole operator 



first. We show results first for Q^ a = 0.3 and the quenching factor Eq. (18) in Fig. |3j 
We can see that the CoGeNT region moves up relative to the DAMA region, so much so 
that it lands above the DAMA region. Agreement can be improved by taking Qno, = 0.45 
as shown in Fig. ([3])b. For the cases that the operator pushes the CoGeNT region above 



the DAMA region, we take Aq = 1 and ro from Eq. (11) since both of these choices tend 



to push CoGeNT down somewhat relative to DAMA (in comparison to Aq = 0.85 and tq 



from Eq. (12)). While the agreement between the two regions is poor, agreement of either 
experiment with the constraints from XENON10 and CDMS-Si is improved, especially for 
DAMA. The results of experiments sensitive to large spin-dependent cross-sections, such as 
COUPP [35] and PICASSO [36], may also be relevant since both the anapole and magnetic 
dipole scattering cross-sections have contributions from a spin- dependent interaction. We 
find, however, that neither of these experiments significantly constrain the otherwise allowed 
parameter space, so we do not show these results on our plots. 

We consider in Fig. Q the effect of the magnetic dipole operator on the scattering regions 
of DAMA and CoGeNT, taking A = 100 MeV. We can see that of all the operators optimal 
agreement between DAMA and CoGeNT is obtained, and improved agreement with the null 
results of XENON10 can be obtained. In addition, we note the effect of changing vq to 270 
km/s is to gain marginal improvement over vq = 220 km/s in terms of the compatibility 
of the regions with each other and with the results of the null experiments as shown in 
Figs. (§,(§. 

Since agreement is optimal between DAMA, CoGeNT and the results of the null 
experiments for the magnetic dipole operator, we show in Fig. ([7]) a sample spectrum 
generated for the magnetic dipole operator with m x = 6.4 GeV, a = 1.9 x 10~ 33 cm 2 . The 
model reproduces very well the DAMA spectrum, but struggles to obtain a large enough 
rise in CoGeNT. The reason for this is the dramatic drop in efficiencies at the low recoil 
energy. The uncorrected rate continues to rise. CoGeNT additionally has data on events in 
one extra bin below their energy threshold, with 26 events, which they discount. This must 
either be explained by a large boundary effect or large errors in the efficiency near threshold. 
Neglecting this, both fits are quite good; Ax 2 = x 2 — Xmm is 2.7 for both CoGeNT and 
DAMA (90% agreement is Ax 2 = 4.91). xLi is 21 - 78 for 17 bins for DAMA and 15.16 for 
27 bins for CoGeNT. It also remains to be explored whether such models could be responsible 
for the 32 events reported by CRESST-II in their oxygen band [37], though it is difficult to 
compare concrete models with no report of the total exposure or efficiency after cuts. 

We have shown here than when experimental uncertainties are taken into account 
appropriately, that a region of parameter space where DAMA and CoGeNT are consistent 
with the results of null experiments is available. Agreement can be improved by considering 
the effects of the scattering primarily through the magnetic dipole operator. We now turn 
to a brief discussion of models where the anapole and magnetic dipole operators are the 
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FIG. 1: left panel: Allowed regions (90 and 99% C.L., corresponding to purple and blue) for 
standard spin-independent scattering, Qjy a = 0.3, Qce from Eq. (17). DAMA regions are shown 



in a darker color than the CoGeNT regions. A green band shows 90% exclusion regions from 
XENON10 depending on the extrapolation of £ e ff below threshold (central values of [15] are taken 
and extrapolated to remain constant (light dashed) below threshold, or to drop linearly to zero 
(dark dashed); these extrapolations correspond roughly to Case 1 and Case 2 of [32]). CDMS-Si 
(red dot-dashed) and SIMPLE (short dashed) constraints are also shown, right panel: Same as left 
panel, but with Qno, = 0.45 and Qce from Eq. (18). 



dominant form of scattering. 



IV. MODELS 



Though different in detail, both the anapole and magnetic dipole operators are velocity 
and momentum suppressed, and thus need sufficiently large cross-sections to explain the 
event rates seen at DAMA and CoGeNT. Therefore the mass of the dark photon A^ that 
mediates the interaction should be fairly light. For example, consider the mass of the 
mediator necessary to generate the large cross-sections for scattering through the anapole 
interaction, Eq. (fTl). That cross-section scales as 
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FIG. 2: left panel: Allowed regions (90 and 99% C.L., corresponding to purple and blue) for standard 



spin-independent scattering, Qno, = 0.45, Qce from Eq. (18) with Aq = 0.85. DAMA regions are 
shown in a darker color than the CoGeNT regions. A green band shows 90% exclusion regions 
from XENON10 depending on the extrapolation of £ e ff below threshold (central values of [15] are 
taken and extrapolated to remain constant (light dashed) below threshold, or to drop linearly to 
zero (dark dashed)). CDMS-Si (red dot-dashed) and SIMPLE (short dashed) constraints are also 
shown, right panel: Same as left panel, but with form factor from Eq. (12) and 20% threshold 
uncertainty in CDMS-Si taken into account. 



where 10~ 34 cm 2 is approximately the size needed for a, from Fig. [3] A cross-section of 
this size is not difficult to generate. Consider a model with a Weyl fermions x an d complex 
scalar <fi with charges +1, —2, respectively. A second Weyl fermion \ c is present for anomaly 
cancellation, but otherwise plays no role; we impose a Z 2 symmetry under which x is °dd 
and all other fields are even in order to restrict the possible interactions of x° with x- The 
allowed renormalizable interactions within the dark sector are then the following: 

£ D X^D, X + X C ^D,X C + \D^\ 2 + V(\4>\ 2 ) + \4>xx + \'4>*X C X° + h.c. (21) 

We assume that <fi dynamically obtains a vev v = {(f)) ~ 10 GeV, which gives mass to x anci 
the dark gauge boson. The x mass term after symmetry-breaking is Majorana, which gives 
a simple explanation for why the anapole operator dominates: the leading vector operator 
Xl^xN^v-N vanishes for Majorana fermions. The dark Majorana particles can couple to the 
dark force because the dark gauge group is broken. Other interactions of x with the dark 
force are higher dimensional and therefore suppressed relative to the anapole interaction. 
Parity is badly broken in the dark sector, and a dark electric dipole moment (EDM), while 
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FIG. 3: left panel: Allowed regions (90 and 99% C.L.) for scattering through the anapole operator, 



QNa = 0.3, Qce from Eq. (17) with Aq = 1, and form factor from Eq. (11). DAMA regions are 
shown in a darker color than the CoGeNT regions. A green band shows 90% exclusion regions 
from XENON10 depending on the extrapolation of £ e g below threshold (central values of [15J are 
taken and extrapolated to remain constant (light dashed) below threshold, or to drop linearly to 
zero (dark dashed)). CDMS-Si (red dot-dashed) and SIMPLE (short dashed) constraints are also 
shown, right panel: Same as left panel, but with QNa = 0.45 and 20% threshold uncertainty in 
CDMS-Si taken into account. 



higher- dimensional, has the same q 2 suppression as the anapole interaction. However, such 
a dark EDM must be generated radiatively and therefore be phase-space suppressed relative 
to the anapole. 

The dark sector then interacts with the Standard Model through kinetic mixing e of the 
light dark force with field strength / M „ with hypercharge: 



CD--f,J^ + ef, u B^. 



^J»UJ 



Then we find 



a 



l(T 34 cm 2 



2 x 10- 



2 ZlOOMeVV 

V m M ) 



8GeVY 



(22) 



(23) 



where m M = \/2g D v is the gauged messenger mass and go is the dark gauge coupling. These 
choices for e and tum are consistent with the bounds on kinetic mixing. 

The magnetic dipole operator may also be easily generated with a sufficiently large cross- 
section. It arises quite naturally when the DM is a Dirac fermion composite. Consider the 
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FIG. 4: left panel: Allowed regions (90 and 99% C.L.) for scattering through the magnetic dipole 
operator, Qn<i = 0.3, Qce from Eq. (17) with An = 1, and form factor from Eq. (11). DAMA 



regions are shown in a darker color than the CoGeNT regions. A green band shows 90% exclusion 
regions from XENON10 depending on the extrapolation of C e s below threshold (central values 
of |15j are taken and extrapolated to remain constant (light dashed) below threshold, or to drop 
linearly to zero (dark dashed)). CDMS-Si (red dot-dashed) and SIMPLE (short dashed) constraints 
are also shown, right panel: Same as left panel, but with Qno, = 0.45 and 20% threshold uncertainty 
in CDMS-Si taken into account. 



case where the dark matter x is a fermionic bound state with compositeness scale A. If x 
is not charged under the dark gauge force, but its constituents are (similar to the neutron 
and the electromagnetic force), then at low momentum transfer its interactions with the 
Standard Model will shut off. We require that the scale A be above the momentum transfer 
~ 30 MeV relevant for scattering at DAMA, so that we may parametrize the interactions 
of x with the dark gauge field in terms of the lowest dimensional gauge-invariant operator, 
Xcr^xF/iv/A. Here we are assuming that parity is not violated in the dark sector, so that 
a dark EDM cannot be generated. In our constraint plots, we have already taken A = 100 
MeV, since this is approximately the minimum that the compositeness scale can be and still 
give a reliable effective theory. With this choice, we find 



a 



4/4 



1.5 x 10 



- 32 cm 2 



nm 



M 



2 x 10" 3 



2 /600MeV\ 4 

V m M J 



(24) 



The additional q 2 dependence in the cross-section that arises from the dark magnetic 
dipole moment tends to lower the direct detection rates at low energies. This allows lighter 
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FIG. 5: left panel: Same as Fig. (pjb (standard WIMP coupling), but with vo = 270 km/s. right 
panel: Same as Fig. {]3J)b (scattering through anapole operator), but with vq = 270 km/s. 
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FIG. 6: Same as Fig. Qb (scattering through magnetic dipole operator), but with vq = 270 km/s. 
One can see that optimal agreement of DAMA and CoGeNT with each other and with the results of 
the null experiments is obtained for scattering through this operator with this set of astrophysical 
and experimental parameters (most importantly, Qno, = 0.45 here). 



DM, since the main obstacle to taking dark matter to be very light is the rapidly rising 
spectrum at low energies in DAMA and CoGeNT, ruining the fit at low recoils. The 
q 2 dependence tends to counteract this, allowing for successful models which evade the 
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FIG. 7: The spectra observed at DAMA (left) and CoGeNT (right) for a benchmark magnetic 
dipole scattering with m x = 6.4 GeV, a = 1.9 x 10 -33 cm 2 . The rate is corrected for CoGeNT 
efficiencies, and these efficiencies are the reason for the drop in the rate in the lowest bin; the 
uncorrected rate continues to rise there. In the right panel, the black solid line is the total predicted 
event rate including signal and background, whereas the red dashed line shows just the signal rate. 



constraints of XENON10, even taking the stochasticity into account. 



V. SUMMARY 

We have studied the implications of experimental uncertainties and non-standard velocity 
dependence on the agreement of CoGeNT and DAMA both with each other and with the 
results of the null experiments. While some marginal agreement can be obtained between 
the two experiments when the sodium quenching factor is pushed to Qn<i = 0.45 and the 
CDMS-Si energy threshold is assumed to have a systematic error of 20% (assumed to be too 
low), optimal achievement between the two regions is not obtained for the standard spin- 
independent case. Agreement can be improved by choosing a different velocity dependence, 
and in particular the magnetic dipole operator gives optimal agreement with all experiments. 
Simple models were constructed where velocity and momentum dependent cross-sections are 
expected to dominate, and it was shown that acceptably large cross-sections can be obtained. 

While the DAMA and CoGeNT signals can be consistent with each other and the 
results of null experiments, large theoretical and experimental uncertainties limit our current 
understanding of the signals and their consistency with models. As we are learning, in an 
age of DM discovery, we must systematically quantify errors on theoretical and experimental 
parameters to determine whether a given DM model is consistent with the signals and with 
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the results of null experiments. Until these uncertainties are reduced and the low mass signals 
can be strongly excluded given the uncertainties, it appears the light DM candidate cannot 
be ruled out. We look forward to further pursuing theoretically well-motivated models of 
low mass DM. 
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Note added: While this work was nearing completion, Refs. [Ml EH] appeared which 
also explore the magnetic dipole operator. Additionally, after this work was submitted to 
arXiv, we became aware of a talk [40J discussing preliminary results from a new method of 
analyzing XENON10 data that would allow it to probe lower recoil energies. Once finalized, 
this method will likely have implications for the scenarios discussed in this paper. 



VI. APPENDIX 

Noting that v -v e = vv e cosx and d 3 v = 27rv 2 dvdcosx, one finds for the standard velocity 
integral in Eq. ^: 



I 



& V Q.gi\y\X esc \%ei-Emin))~T~ ^j-PSI^y^esc •Emini^ei •Eesci%min~r-Ee)y 



\v\>v min v N N 



where 



(25) 



v.2 / _ o X~ 



asi = Vn (erf(a; m i n + x e ) - erf(x min - x e )) - 4x e e x%sc I 1 + x 2 esc - -^ - a^J j ,(26) 
Psi = V^ (erf (aw) - eii(x min - x e )) (27) 



with 

A . , "1 ,- „, _„2 /2x 3 

JM = Ax e v 



y/ir erf(x esc ) - e x ^x esc I — ^ + 1 



(29) 



Likewise one can compute the relevant integral in Eq. (p3J) for v 2 dependent cross-sections. 
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It is 

" J \yiV e )(l V ——0L V 2 \J\X esc \X e -rX m i n ) )~r ——p v 2KyyX esc X m i n -rX e , X esc' X min* X e ) j 

W\>v mm ^ ^ 

(30) 
where 



'«,.j --e " [4^n(Z.r; + [)VYn,r l , ini .r, l )e^ min+Xe)2 



-3^ (2Xg + l) erf(x min + x e )e (Xmm+Xe)2 - Qx e e iXminXe + 24Xg SC x e e" x - c+(Xmsil+:re)2 

-PiT ■ p^ X min x e _l_ f?™ . f?™ _l_ Q T -Jp- ;E esc + (:C m .i n +£ e J _i_ OJ™ f D~ x esc'\ x 'min~r x e) 



— e esc ( — H - 2x e (x min — x esc j I (31) 



/3„2 = 2e Xesc 



o ( x Ln - Ow + O 3 ) ~ ^ e *""" "'"^^ (v^ ( 2 ^e + l ) e x ™ in+xl (erf(x min - x e 



K:LL \ J ^escj I z, ° ^min ^ "^e/ ' ^y^esc ' ^-^ej^ J q^. ^ A e T iu *e esc 

+ 10x e (^£ esc + x min j + 15x e (x esc — x min j + <lx esc — L\jx esc x min + 0£ m j n JJ [oZ) 



The standard velocity integral Eq. (25) (weighted with q 2 or q 4 ) can be combined with 
the v 2 integral to give the total rate for the scattering through the anapole and magnetic 
dipole operators. 
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